El modelado matemático en el estudio de mecanismos de señalización celular: apoptosis inducida por interferón de tipo I

Comunicaciones a congreso

Accurate cell segmentation in microscopy images using membrane patterns

Motivation: Identifying cells in an image (cell segmentation) is essential for quantitative single-cell biology via optical microscopy. Although a plethora of segmentation methods exists, accurate segmentation is challenging and usually requires problem-specific tailoring of algorithms. In addition, most current segmentation algorithms rely on a few basic approaches that use the gradient field of the image to detect cell boundaries. However, many microscopy protocols can generate images with characteristic intensity profiles at the cell membrane. This has not yet been algorithmically exploited to establish more general segmentation methods. Results: We present an automatic cell segmentation method that decodes the information across the cell membrane and guarantees optimal detection of the cell boundaries on a per-cell basis. Graph cuts account for the information of the cell boundaries through directional cross-correlations, and they automatically incorporate spatial constraints. The method accurately segments images of various cell types grown in dense cultures that are acquired with different microscopy techniques. In quantitative benchmarks and comparisons with established methods on synthetic and real images, we demonstrate significantly improved segmentation performance despite cell-shape irregularity, cell-to-cell variability and image noise. As a proof of concept, we monitor the internalization of green fluorescent protein-tagged plasma membrane transporters in single yeast cells. Availability and implementation: Matlab code and examples are available at http://www.csb.ethz.ch/tools/cellSegmPackage.zip. Contact: [email protected] or [email protected] Supplementary information: Supplementary data are available at Bioinformatics onlin

MetaNetX.org: a website and repository for accessing, analysing and manipulating metabolic networks

Summary: MetaNetX.org is a website for accessing, analysing and manipulating genome-scale metabolic networks (GSMs) as well as biochemical pathways. It consistently integrates data from various public resources and makes the data accessible in a standardized format using a common namespace. Currently, it provides access to hundreds of GSMs and pathways that can be interactively compared (two or more), analysed (e.g. detection of dead-end metabolites and reactions, flux balance analysis or simulation of reaction and gene knockouts), manipulated and exported. Users can also upload their own metabolic models, choose to automatically map them into the common namespace and subsequently make use of the website's functionality. Availability and implementation: MetaNetX.org is available at http://metanetx.org. Contact: [email protected]

Multistable and dynamic CRISPRi-based synthetic circuits

Gene expression control based on CRISPRi (clustered regularly interspaced short palindromic repeats interference) has emerged as a powerful tool for creating synthetic gene circuits, both in prokaryotes and in eukaryotes; yet, its lack of cooperativity has been pointed out as a potential obstacle for dynamic or multistable synthetic circuit construction. Here we use CRISPRi to build a synthetic oscillator ("CRISPRlator"), bistable network (toggle switch) and stripe pattern-forming incoherent feed-forward loop (IFFL). Our circuit designs, conceived to feature high predictability and orthogonality, as well as low metabolic burden and context-dependency, allow us to achieve robust circuit behaviors in Escherichia coli populations. Mathematical modeling suggests that unspecific binding in CRISPRi is essential to establish multistability. Our work demonstrates the wide applicability of CRISPRi in synthetic circuits and paves the way for future efforts towards engineering more complex synthetic networks, boosted by the advantages of CRISPR technology

Efficient characterization of high-dimensional parameter spaces for systems biology

BACKGROUND: A biological system's robustness to mutations and its evolution are influenced by the structure of its viable space, the region of its space of biochemical parameters where it can exert its function. In systems with a large number of biochemical parameters, viable regions with potentially complex geometries fill a tiny fraction of the whole parameter space. This hampers explorations of the viable space based on "brute force" or Gaussian sampling. RESULTS: We here propose a novel algorithm to characterize viable spaces efficiently. The algorithm combines global and local explorations of a parameter space. The global exploration involves an out-of-equilibrium adaptive Metropolis Monte Carlo method aimed at identifying poorly connected viable regions. The local exploration then samples these regions in detail by a method we call multiple ellipsoid-based sampling. Our algorithm explores efficiently nonconvex and poorly connected viable regions of different test-problems. Most importantly, its computational effort scales linearly with the number of dimensions, in contrast to "brute force" sampling that shows an exponential dependence on the number of dimensions. We also apply this algorithm to a simplified model of a biochemical oscillator with positive and negative feedback loops. A detailed characterization of the model's viable space captures well known structural properties of circadian oscillators. Concretely, we find that model topologies with an essential negative feedback loop and a nonessential positive feedback loop provide the most robust fixed period oscillations. Moreover, the connectedness of the model's viable space suggests that biochemical oscillators with varying topologies can evolve from one another. CONCLUSIONS: Our algorithm permits an efficient analysis of high-dimensional, nonconvex, and poorly connected viable spaces characteristic of complex biological circuitry. It allows a systematic use of robustness as a tool for model discrimination

Topological augmentation to infer hidden processes in biological systems

Motivation: A common problem in understanding a biochemical system is to infer its correct structure or topology. This topology consists of all relevant state variables—usually molecules and their interactions. Here we present a method called topological augmentation to infer this structure in a statistically rigorous and systematic way from prior knowledge and experimental data. Results: Topological augmentation starts from a simple model that is unable to explain the experimental data and augments its topology by adding new terms that capture the experimental behavior. This process is guided by representing the uncertainty in the model topology through stochastic differential equations whose trajectories contain information about missing model parts. We first apply this semiautomatic procedure to a pharmacokinetic model. This example illustrates that a global sampling of the parameter space is critical for inferring a correct model structure. We also use our method to improve our understanding of glutamine transport in yeast. This analysis shows that transport dynamics is determined by glutamine permeases with two different kinds of kinetics. Topological augmentation can not only be applied to biochemical systems, but also to any system that can be described by ordinary differential equations. Availability and implementation: Matlab code and examples are available at: http://www.csb.ethz.ch/tools/index. Contact: [email protected]; [email protected] Supplementary information: Supplementary data are available at Bioinformatics onlin

Effect of Immunosuppression on T-Helper 2 and B-Cell Responses to Influenza Vaccination

Background. Influenza vaccine immunogenicity is suboptimal in immunocompromised patients. However, there are limited data on the interplay of T- and B- cell responses to vaccination with simultaneous immunosuppression. Methods. We collected peripheral blood mononuclear cells from transplant recipients before and 1 month after seasonal influenza vaccination. Before and after vaccination, H1N1-specific T- and B-cell activation were quantified with flow cytometry. We also developed a mathematical model using T- and B-cell markers and mycophenolate mofetil (MMF) dosage. Results. In the 47 patients analyzed, seroconversion to H1N1 antigen was demonstrated in 34%. H1N1-specific interleukin 4 (IL-4)-producing CD4+ T-cell frequencies increased significantly after vaccination in 53% of patients. Prevaccine expression of H1N1-induced HLA-DR and CD86 on B cells was high in patients who seroconverted. Seroconversion against H1N1 was strongly associated with HLA-DR expression on B cells, which was dependent on the increase between prevaccine and postvaccine H1N1-specific IL-4+CD4+ T cells (R2 = 0.35). High doses of MMF (≥2 g/d) led to lower seroconversion rates, smaller increase in H1N1-specific IL-4+CD4+ T cells, and reduced HLA-DR expression on B cells. The mathematical model incorporating a MMF-inhibited positive feedback loop between H1N1-specific IL-4+CD4+ T cells and HLA-DR expression on B cells captured seroconversion with high specificity. Conclusions. Seroconversion is associated with influenza-specific T-helper 2 and B-cell activation and seems to be modulated by MM

A synthetic low-frequency mammalian oscillator

Circadian clocks have long been known to be essential for the maintenance of physiological and behavioral processes in a variety of organisms ranging from plants to humans. Dysfunctions that subvert gene expression of oscillatory circadian-clock components may result in severe pathologies, including tumors and metabolic disorders. While the underlying molecular mechanisms and dynamics of complex gene behavior are not fully understood, synthetic approaches have provided substantial insight into the operation of complex control circuits, including that of oscillatory networks. Using iterative cycles of mathematical model-guided design and experimental analyses, we have developed a novel low-frequency mammalian oscillator. It incorporates intronically encoded siRNA-based silencing of the tetracycline-dependent transactivator to enable the autonomous and robust expression of a fluorescent transgene with periods of 26 h, a circadian clock-like oscillatory behavior. Using fluorescence-based time-lapse microscopy of engineered CHO-K1 cells, we profiled expression dynamics of a destabilized yellow fluorescent protein variant in single cells and real time. The novel oscillator design may enable further insights into the system dynamics of natural periodic processes as well as into siRNA-mediated transcription silencing. It may foster advances in design, analysis and application of complex synthetic systems in future gene therapy initiatives

BaSysBio: Towards an understanding of dynamic transcriptional regulation at global scale in bacteria: a systems biology approach

9th International Conference on Systems Biology, ICSB 2008, 22-28 August, Goteborg, SwedenBaSysBio adopts a systems biology approach in which quantitative experimental data will be generated for each step of the information flow and will fuel computational modelling. High throughput technologies (living cell arrays, tiling DNA microarrays, MDLC proteomics and quantitative metabolomics) will be developed in conjunction with new computational modelling concepts to facilitate the understanding of biological complexity. Models will simulate the cellular transcriptional responses to environmental changes and their impact on metabolism and proteome dynamics. The iterative process of simulations and model-driven targeted experiments will generate novel hypotheses about the mechanistic nature of dynamic cellular responses, unravel emerging systems properties and ultimately provide an efficient roadmap to tackle novel, pathogenic organismsN

SBML Level 3: an extensible format for the exchange and reuse of biological models

Abstract Systems biology has experienced dramatic growth in the number, size, and complexity of computational models. To reproduce simulation results and reuse models, researchers must exchange unambiguous model descriptions. We review the latest edition of the Systems Biology Markup Language (SBML), a format designed for this purpose. A community of modelers and software authors developed SBML Level 3 over the past decade. Its modular form consists of a core suited to representing reaction‐based models and packages that extend the core with features suited to other model types including constraint‐based models, reaction‐diffusion models, logical network models, and rule‐based models. The format leverages two decades of SBML and a rich software ecosystem that transformed how systems biologists build and interact with models. More recently, the rise of multiscale models of whole cells and organs, and new data sources such as single‐cell measurements and live imaging, has precipitated new ways of integrating data with models. We provide our perspectives on the challenges presented by these developments and how SBML Level 3 provides the foundation needed to support this evolution